WHAT IS CLAIMED IS: 

1 . A method of modeling flame propagation comprising: 

defining a flame surface area density of a flame as a flame surface area per unit 
volume of the flame; 

expressing flame progress as generation of the flame surface area density in terms 
of at least one of a turbulent combustion and a laminar combustion; 

determining flame growth resulting from turbulent combustion as being inversely 
proportional to a chemical reaction characteristic time and as a function of a turbulent 
Reynolds number; and 

modeling the flame propagation based on the flame growth. 

2. The flame propagation modeling method as recited in claim 1 , further 
comprising 

determining the flame growth resulting from laminar combustion as being 
proportional to both a laminar flame speed and to a ratio of a temperature of a burned 
portion to a temperature of an unburned portion and as a function of the Karlowitz number. 

3. The flame propagation modeling method as recited in claim 1, wherein 
the generation of the flame surface area density is expressed as a combination of 

the turbulent combustion and the laminar combustion. 

4. The flame propagation modeling method as recited in claim 1 , wherein 
the flame growth resulting from the turbulent combustion is calculated based on 

the flame growth being inversely proportional to the chemical reaction characteristic time 
and proportional to both the turbulent Reynolds number raised to an exponential power 
and a stretch rate of the flame. 

5. The flame propagation modeling method as recited in claim 2, wherein 
the flame generation is further expressed as transport of the flame surface area 

density, which is expressed in terms of flame growth resulting from turbulent combustion 
and flame growth resulting from laminar combustion; and 
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the flame growth resulting from laminar combustion being expressed as 
proportional to the laminar flame speed, to the ratio of the temperature of a burned portion 
to the temperature of an unbumed portion, and to an exponential function of the Karlowitz 
number. 

5 

6. The flame propagation modeling method as recited in claim 5, wherein 
the exponential function of the Karlowitz number is the base of the natural 
logarithm raised to the power of the Karlowitz number. 

10 7. The flame propagation modeling method as recited in claim 1 , wherein 

the flame growth resulting from the turbulent combustion is expressed as follows: 

S T =aj(Re t ) a2 r- E, 

K 

where St represents flame growth resulting from turbulent combustion, 
E represents flame surface area density, k represents turbulence strength, 
15 e represents turbulence dissipation rate, Re t represents turbulent Reynolds number, r 
represents flame stretch rate, and a } and a 2 are model constants. 

8. The flame propagation modeling method as recited in claim 2, wherein 
the flame growth resulting from the laminar combustion is expressed as follows: 

20 S L = fit exp( - p 2 Ka) — U L I 2 , 

Tu 

where S L flame growth resulting from laminar combustion, 
E represents flame surface area density, U L represents laminar flame speed, 
T b represents burned gas temperature, T u represents unburned gas temperature, Ka 
represents Karlowitz number, and /?; and J3 2 are model constants. 

25 

9. The flame propagation modeling method as recited in claim 1 , wherein 
the flame generation is further expressed as transport of the flame surface area 

density, which is expressed in terms of flame growth resulting from turbulent combustion 
and flame growth resulting from laminar combustion; and 
30 the flame generation is suppressed by a resistance force imposed by air. 
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10. The flame propagation modeling method as recited in claim 1 , wherein 
transport, generation, and diffusion of the flame surface area density are expressed 
as follows: 



9Z du ,Z 3 



+ ai (Ref 2 r-Z + fa exp(-frKa) —U L E 2 - D, 



^C7 C dx i J K Tu 



5 where 27 represents flame surface area density, k represents turbulence strength, e 

represents turbulence dissipation rate, Re t represents turbulent Reynolds number, r 
represents flame stretch rate, U L represents laminar flame speed, 7& represents burned gas 
temperature, T u represents unburned gas temperature, Ka represents Karlowitz number, v, 
represents turbulent kinematic viscosity, a c represents turbulent Schmidt number, 
10 D represents air resistance force, and a/, <*2, fti and are model constants. 

11. A method of modeling flame propagation comprising: 

defining a flame surface area density of a flame as a flame surface area per unit 
volume of the flame; 

15 expressing flame progress as generation of the flame surface area density in terms 

of at least one of a turbulent combustion and a laminar combustion; 

determining flame growth resulting from laminar combustion as being proportional 
to both a laminar flame speed and to a ratio of a temperature of a burned portion to a 
temperature of an unburned portion and as a function of the Karlowitz number; and 

20 modeling the flame propagation based on the flame growth. 

12. The flame propagation modeling method as recited in claim 11, wherein 
the generation of the flame surface area density is expressed as a combination of 

the turbulent combustion and the laminar combustion. 

25 

13. The flame propagation modeling method as recited in claim 12, further 
comprising 

determining flame growth resulting from the turbulent combustion is calculated 
based on a flame growth being inversely proportional to a chemical reaction characteristic 
30 time and proportional to both a turbulent Reynolds number raised to an exponential power 
and a stretch rate of the flame. 
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14. The flame propagation modeling method as recited in claim 11, wherein 
the flame generation is further expressed as transport of the flame surface area 

density, which is expressed in terms of flame growth resulting from turbulent combustion 
5 and flame growth resulting from laminar combustion; and 

the flame growth resulting from laminar combustion being expressed as 
proportional to the laminar flame speed, to the ratio of the temperature of a burned portion 
to the temperature of an unbumed portion, and to an exponential function of the Karlowitz 
number. 

0 

15. The flame propagation modeling method as recited in claim 14, wherein 
the exponential function of the Karlowitz number is the base of the natural 

logarithm raised to the power of the Karlowitz number. 



15 16. The flame propagation modeling method as recited in claim 13, wherein 

the flame growth resulting from the turbulent combustion is expressed as follows: 

ST^oaiRet)* 2 T- 27, 

K 

where St represents flame growth resulting from turbulent combustion, 
27 represents flame surface area density, k represents turbulence strength, 
20 e represents turbulence dissipation rate, Re t represents turbulent Reynolds number, r 
represents flame stretch rate, and oa and a 2 are model constants. 



17. The flame propagation modeling method as recited in claim 16, wherein 
the flame growth resulting from the laminar combustion is expressed as follows: 

25 S L = fij exp( - feKa) — U L 27 2 , 

Tu 

where Sl flame growth resulting from laminar combustion, 
27 represents flame surface area density, Ul represents laminar flame speed, 
T b represents burned gas temperature, T u represents unbumed gas temperature, Ka 
represents Karlowitz number, and /?; and y8 2 are model constants. 

30 
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18. The flame propagation modeling method as recited in claim 11, wherein 
the flame growth resulting from the laminar combustion is expressed as follows: 

S L = faexp(-fi 2 Ka)'^ U L Z 2 f 

Tu 

where Sl flame growth resulting from laminar combustion, 
5 Z represents flame surface area density, Ul represents laminar flame speed, 

Tb represents burned gas temperature, T u represents unbumed gas temperature, Ka 
represents Karlowitz number, and /?/ and f} 2 are model constants. 

19. The flame propagation modeling method as recited in claim 1 1 , wherein 
10 the flame generation is further expressed as transport of the flame surface area 

density, which is expressed in terms of flame growth resulting from turbulent combustion 
and flame growth resulting from laminar combustion; and 

the flame generation is suppressed by a resistance force imposed by air. 

15 20. The flame propagation modeling method as recited in claim 1 1 , wherein 

transport, generation, and diffusion of the flame surface area density are expressed 
as follows: 



dt dx, dx ; 



{<T C dx r 



+ a } (Ref 2 T-Z+fa exp(-fi 2 Ka) — U L Z 2 - D, 

K Tu 



where Z represents flame surface area density, k represents turbulence strength, e 
20 represents turbulence dissipation rate, Re t represents turbulent Reynolds number, r 

represents flame stretch rate, U L represents laminar flame speed, T b represents burned gas 
temperature, T u represents unbumed gas temperature, Ka represents Karlowitz number, v t 
represents turbulent kinematic viscosity, a c represents turbulent Schmidt number, 
D represents air resistance force, and ay, a 2f fii and fi 2 are model constants. 

25 



- 16- 



